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ABSTRACT 

In this work we use a sample of 318 radio-quiet quasars (RQQ) to investigate the dependence of the 
ratio of optical/UV flux to X-ray flux, a ox , and the X-ray photon index, Tx, on black hole mass, UV 
luminosity relative to Eddington, and X-ray luminosity relative to Eddington. Our sample is drawn 
from the literature, with X-ray data from ROSAT and Chandra, and optical data mostly from the 
SDSS; 153 of these sources have estimates of Tx from Chandra. We estimate Mbh using standard 
estimates derived from the H/3, Mg II, and C IV broad emission lines. Our sample spans a broad 
range in black hole mass (10 6 < M B h/M q < 10 10 ), redshift (0 < z < 4.8), and luminosity (10 43 < 
ALA(2500A)[erg s" 1 ] < 10 48 ). We find that a ox increases with increasing Mbh and Luv / LEdd, and 
decreases with increasing Lx/LEdd- In addition, we confirm the correlation seen in previous studies 
between Tx and Mbh and both Luv I LEdd and Lx/LEdd] however, we also find evidence that the 
dependence of Tx of these quantities is not monotonic, changing sign at Mbh ~3x 1O 8 M . We argue 
that the a ox correlations imply that the fraction of bolometric luminosity emitted by the accretion 
disk, as compared to the corona, increases with increasing accretion rate relative to Eddington, rh. 
In addition, we argue that the Tx trends are caused by a dependence of X-ray spectral index on rh. 
We discuss our results within the context of accretion models with comptonizing corona, and discuss 
the implications of the a ox correlations for quasar feedback. To date, this is the largest study of the 
dependence of RQQ X-ray parameters on black hole mass and related quantities, and the first to 
attempt to correct for the large statistical uncertainty in the broad line mass estimates. 
Subject headings: accretion disks — galaxies: active — quasars: general — ultraviolet: galaxies — 
X-rays: galaxies 



1. INTRODUCTION 

The extraordinary activity associated with quasars 
involves accretion onto a supermassive black hole 
(SMBH), with the UV/optical emission arising from 
a ge ometrically thin, optically thick cold accretion 
disk (Sha kura fc Svunvaevi Il973h . and the X-ray con- 
tinuum arising from a hot, optically thin corona 
that Compton upscatters the disk UV photons (e.g., 
lHaardt k Maraschil 1 1 9 9 lh . In highly accreting ob- 
jects, like quasars (0.01 < L^JLtem "> 1» e -K- 



Woo fc Urr^l2002tlVestergaaTdll2004lMcLTire fc Dunlod 



20041: iKollmeier et al l l2006t) . the X-ray plasma ge- 
ometry is expected to be that of a hot, pos- 
sibly patchy, io nized 'skin' that sandwiches the 
cold disk (e.g., IBisnovatvi-Kogan fc Blinnikovl 119771 : 



iLiang fc Pricelll977l : lNavakshinll2000f) . However, the ev- 
idence for this is not conclusive, and relies on data from 
X-ra y binaries and low- z sources (e.g., see the dicussion 
by ICzernv et al.l I2003D . Other geometries are possi- 
ble, including an ac cretion disk that evaporates into a 
hot i nner flow (e.g., IShapiro et al.l 119761 : IZdziarski et al.l 
fl999T) . or a combination of a hot inner flow and a corona 
that sandwiches the di sk (e.g., iPoutanen et al.l [19971 : 
ISobolewska et al.l l2004af l. Furthermore, radiation pres- 

Electronic address: bkelly@as.arizona.edu 

1 Dept. of Physics and Astronomy, Robinson Hall, Tufts Uni- 
versity, Medford, MA 02155 (present address) 



sure can drive an outflow from the disk into the corona 
if the t wo are cospat ial, thus altering the physics of the 
corona (|Progall2005f ). Investigations of how quasar X- 
ray parameters depend on black hole mass, Mbh, and 
accretion rate relative to Eddington, m, offer important 
constraints on models of the disk/corona system. 

There have been attempts to link the evolution 
of SMBHs to analytic and semi-analytic models of 
structure formation (e.g.. |Kauffmann fc HaehneltJ 120001 : 
lHatziminaoglou et al.l 120031 : iBromlev et al.ll2004T) 7 where 
black holes grow by accreting gas funneled towards 
the center during a galaxy merger until feedback en- 
ergy from the SMBH expels gas and shuts off the ac- 
cretion process (e.g., [Silk & Rccs 1998; Fabian 119991 : 
iWvithe fc Loebl MMi iBegelman k Nathl I2005T 1. This 
'self-regulated' growth of black holes has recently been 
successfully applied in smoothed particle hydrodynamic s 
simulations ()Di Matteo et al.ll2005l : ISpringel et al.ll2005t ). 
Within this framework, the AGN or quasar phase oc- 
curs during the episode of significant accretion that fol- 
lows the galaxy merger, persisting until feedback from 
the b lack hole 'blows' the gas away (e.g.. iHopkins et"aTI 
2006). Hydrodynamic calculations have shown that line 
pressure is more efficient than thermal pressure at driving 
an outflow (Proga 2007), and therefore, the efficiency of 
AGN feedback depends on the fraction of energy emitted 
through the UV/disk component as compared to the X- 
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ray/corona component. If the fraction of energy emitted 
in the UV as compared to the X-ray depends on Mbh 
or m, then it follows that the efficiency of AGN feedback 
will also depend on Mbh and m. This has important 
consequences for models of SMBH growth, as the SMBH 
may become more or less efficient at driving an outflow 
depending on its mass and accretion rate. Studies of the 
dependence of quasar X-ray/UV emission on black hole 
mass and accretion rate are therefore important as they 
allow us to constrain a Mbh- or rri-dependent feedback 
efficiency. 

Numerous previous studies have searched 
for a luminosity and redshift dependence of 
a ox — —0.384 log Lx/Lu vj the ratio of X -r ay to 
UV/optical flux ( e.g.. lAvni fc Tananbauml Il982t 
Wilkes et all 11994 lYuan et all Il998t IVignali et all 



2003b; 
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20051: iSteffen et all 



Strateva et al.l 

20071) . and fx, t h e X-ray spectra slope 



e.g., iReeves fe Turner! 120001: [Bechtold et all 120031: 



Dai et all 12004 iRisaliti fc Elvisl 120051 : iGrupe et all 
20061) . Most studies have found a correlation between 



a ox and UV luminosity, Luv, while the existence of 
a correlation between q nx and z is still a matter of 
debate (e.g. iBechtold et all 120031: IVignali et alj |2003bt 



ISteffen et allboofl: Uust et alJ l200ftlKellv et alJ I2007F . 

In addition, studies of Fx have produced mixed results. 
Some authors have claimed a correlation between fy 



and luminosity (e.g., IBechtold et alj 120031: iDai et al 
20041) or redshift fe-g., [Reeves et alJll997t IVignali et al 



1999b iPage et ah! 12003). while, others find no evidence 
for a correlation between and L\jy or z (e.g. 



Vignali et al.l 12001 IRisaliti k Elvisl 12003: iKellv et all 
2007). 

A correlation between Tx and the FWHM of 
the H/3 line has also b een found (e.g., iBoller et al.l 
119961 : iBrandt et al.lll997l ). suggesting a correlation be- 
tween r^_and black hole mass or Ed dingto n ratio 
(e.g., lLaor et all 119971: IBrandt fe Boiler] I1998J ). Re- 
cently, it has become possible to obtain estimates 
of Mbh for broad line AGN by calibrating re- 
sults from r e verbe ration mapping (jPeterson et alJl2004t 
IKaspi et al.l 120051) for use on single-epoch spectra 
(IWandel et al.lll999l:IVestergaard 2002; M cLure fe Jarvisj 
'2002 t IVestergaard fc Petersonl 120061 : IKellv fc Bechtoldl 
2007|) . This has enabled some authors to confirm a cor- 
relation betwee n Tx and either Mbh or Ltoi/LEdd (e.g.. 



Lu fc Yd Il99l [Gje iiihski & Done 2004; Porauc t et all 



20041: iPiconcelli et a l. 2005; Shc mmer et~aTll2006h . where 



the X-ray continuum hardens with increasing Mbh or 
softens with increasing L^oi/ L E dd- In addition, previ- 
ous work has also found evidence for quasars becom- 
ing more X-ray quiet as Mbh or L b i/L E dd increase 
(jBrunner et al.ll997tlWang et al.ll2004D : however, studies 
involving the dependence of a ox on Mbh ° r Lb i/L E dd 
have remained rare compared to studies of Tx ■ It is im- 
portant to note that the correlations inferred in previous 
work generally employ broad line mass estimates in com- 
bination with a constant bolometric correction. There- 
fore, most of the correlations found in previous work are, 
strictly speaking, between Tx or a ox and the estimates 
M BH oc LlFWHM 2 and L bo i/L Edd oc L 1 ^ 1 FWHM~ 2 , 
where 7 ~ 0.5. 
In this work, we investigate the dependence of a ox and 



Tx on black hole mass, optical/UV luminosity relative 
to Eddington, and X-ray luminosity relative to Edding- 
ton. We combin e the main Sloan Di gital Sky Survey 
(S PSS) sample oflStrateva et al.l (|2005l ) with the sample 
of IKellv et ail (|2007D . creating a sample of 318 radio- 
quiet quasars (RQQ) with X-ray data from ROSAT and 
Chandra, and optical spectra mostly from the SDSS; 
153 of these sources have estimates of Tx from Chan- 
dra. Because the X-ray emission in radio-loud sources 
can have an additional component from the jet (e.g., 
IZamorani et~aIlll981l : IWilkes fc Elvislll987f ). we focus our 
analysis on the radio-quiet majority. Our sample has a 
detection fraction of 87% and spans a broad range in 
black hole mass (10 6 < M BH /M Q < 10 10 ), redshift (0 < 
z < 4.8), and luminosity (10 43 < AL A (2500A)[erg s" 1 ] < 
10 48 ), enabling us to effectively look for trends regarding 
a ox and T X - 

The outline of this paper is as follows. In § [2] wc de- 
scribe the construction of our sample, and in § [3] we 
describe the procedure we used to fit the optical con- 
tinuum and emission lines. In § 0] we describe how we 
obtain broad line mass estimates, our bolometric cor- 
rection, and argue that a constant bolometric correction 
provides a poor estimate of the bolometric luminosity. In 
§ [5] we describe the results from a regression analysis of 
a ox on MBH,Luv/LEdd, and Lx/Lsdd, and in §[6] we 
report evidence for a non-monotonic dependence of Tx 
on either M B h , L uv /L Edd , and L x /L Ed d- In §[7] we dis- 
cuss our results within the context of AGN disk/corona 
models, and we discuss the implications for a dependence 
of quasar feedback efficiency on black hole mass or ac- 
cretion rate. In §[8] we summarize our main results. 

We adopt a cosmology based on the the WMAP 
best-fit parameters (h — 0.71, f2 TO = 0.27, f^A = 0.73, 
ISpergel et all l2003f) . For ease of notation, we de- 
fine L uv = i/i y (2500A),Lx = vL v {2 keV),fov = 
log^i^(2500A), and m B H = log M BH /M Q . 

2. SAMPLE CONSTRUCTION 



In t his analysis we combine 169 RQQs from lKellv et al.l 
(|2007l hereafter K07) with 149 RQQs from the main 
SDSS sample of IStrateva et all (|2005l hereafter S05) 
to create a sample of 318 RQQs. Out of these 318 
sources, 276 (86.8%) are detected in the X-rays. The 
z < 4 sources from the K07 sample were selected 
by cross-correlating t he SDSS DR3 quasar catalogue 
(|Schneider et al.ll2005l ) with the Chandra public archive 
as of 2005 February 22. The z > 4 sources from the 
K07 sample consist of targeted Chandra RQQs taken 
from the lite rature (|Bechtold et al.l l2003t IVignali et al.l 
I200H l2003al) . and new observations reported by K07. 
The sources taken from S05 were selected from the SDSS 
to be contained within the inner 19' of ROSAT PSPC 
pointings with exposure times > 11 ksec. The X-ray data 
for both samples are as reported by S05 and K07. 

Both the S05 and K07 samples consist only of radio- 
quiet quasars. We focus our analysis on the radio- 
quiet majority because the radio-loud sources have an 
additional c omponent of X-ray emission arising from 
the jet (e.g.. IZamorani et al.lll981b IWilkes fc Ehdlll987l : 
iWorrall et alj 119871 ). IrTaddition. both S05 and K07 
omitted BAL QSOs when possible. It is necessary to 
remove the BAL QSOs because their high column den- 
sity gives them the appearance of being X-ray weak (e.g., 
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iGreen et al.l I200H : iGallaeher etlll [200l 12001 . poten- 
tially biasing our analysis. However, neither S05 nor K07 
were able to remove the high-ionization BAL quasars for 
z < 1.5, as their identification requires observations of 
the C IV line. In addition, low-ionization BALs can be 
identified at 0.45 < z < 2.25 based on Mg II absorption. 
iReichard et al.l (|2003[ ) found the fraction of BALs in the 
SDSS to be ~ 14%, and therefore we expect there to be 
25 ± 5 BALs in our sample at z < 1.5. This number may 
be higher if one rel axes the definition of a BAL quasar 
([Trump et al.ll2006l ). 

We exclude five sources from the S05 sample due to 
significant intrinsic narrow UV absorption or obvious 
host galaxy contamination: Source SDSS J103747.4- 
001643.9 [z = 1.500) had significant C IV absorption, 
and sources SDSS J124520.7-002128.1 (z = 2.354) and 
SDSS J103709.8+000235.2 (z = 2.679) had significant 
absorption in both C IV and La. These three sources 
were omitted because the absorption prohibits obtain- 
ing an accurate line width measurement, necessary for 
broad line mass estimates, and to ensure that the X- 
ray emission under study is not effected by the absorp- 
tion. Sources SDSS J230440.6-082220.8 (z = 0.201) and 
SDSS J023306.0+003856.4 (z = 0.244) have a significant 
host-galaxy component in their spectra. In addition, we 
exclude source SDSS J144340.8+585653.2 {z = 4.278) 
from the K07 sample because it has significant UV ab- 
sorption. We removed source SDSS J142414.1+421400.1 
(z = 1.608) from the K07 sample and source SDSS 
J170441.4+604430.5 (PG 1704+608, z = 0.372) from the 
S05 sample, as both sources are radio-loud. 

We could not estimate black hole masses for 
sources SDSS J083206.0+524359.3 (z = 1.573), SDSS 
J144231. 7+011055. 3 (z = 4.507), and PC 0910+5625 
(z = 4.035). All three of these sources are from the 
K07 sample. The region containing the Mg II emission 
line for SDSS J083206.0+524359.3 was missing from the 
SDSS spectrum, the emission lines are too weak for SDSS 
J144231. 7+011055. 3, and an optical spectrum was not 
available for PC 0910+5625. 

3. OPTICAL/UV SPECTRAL FITS 

Optical spectra were obtained for most sources from 
the SDSS. We also obtained spectra for s o me o f 
the high redshift q uasars from | Anderson et al.l d200lT) , 
iPeroux et all (|2001l ). and lConstantin et all (|2002t) . The 
values of Luy and a, L v oc v~ a , for the K07 sources 
are taken from K07. We processed the optical spectra 
for the S05 sources in the same manner as for the K07 
sources. We do this for consistency and because S05 did 
not correct for quasar iron emission. 

3.1. Continuum Fitting 

As described by K07, we corrected the optical spec- 
tra for Galactic absorption u si ng th e E(B — V) val- 
ues taken from ISchlegel et al.l (|1998h . as listed in the 
NASA/IPAC Ex tragalactic Database (NED), and the ex- 
tinction curve of ICardelli et al.l (|1989l ). assuming a value 
of Ay/E(B — V) = 3.1. We model the continuum as 
a power law of the form /„ oc ^~ Q , and the Fe emis- 
sion as a scaled and broadened iron template extracted 
fro m I Zw I. The op t ical i ron template was extracted 
bv lVeron-Cettv et all (120041). and the U V iron template 
was extracted bv lVestergaard fe Wilkesl (|2001l ). The con- 



tinuum and iron emission were fit simultaneously us- 
ing the Levenberg-Marquardt method for nonlinear \ 2 ~ 
minimization. Continuum flux densities were then esti- 
mated using the power law parameters. 

We were not able to use a power-law fit to calculate 
Ljjy for the z < 0.4 sources, as the SDSS spectral range 
for these sources does not contain the rest-frame UV con- 
tinuum. Instead, we use the luminosity of the broad 
component of the H/3 emission line, Lh/3, as a proxy for 
Ljjy. It is preferable to use the broad H/3 emission line 
luminosity over, say, the optical continuum luminosity as 
a proxy for Ljjy because the broad H/3 emission line is 
not contaminated by emission from the host galaxy, and 
thus should provide an approximately unbiased estimate 
of Luy. Host-galaxy contamination is likely negligible 
for all z > 0. 4 sources, as vL* ~ 1 44 ergs s -1 at 2500A 
for galaxies (|Budavari et al.l [20051 ). We fit a power-law 
relationship between Luy and Lh/3 using the 44 sources 
at 0.4 < z < 0.9 for which a measurement of both quan- 
tit ies if availab le. We used the linear regression method 
of iKellvl (|2007l) . which allows for measurement errors in 
both variables, and find 



AL A (2500A) 
10 44 ergs s" 1 



(1.556+0.282) 



L 



H/3 



10 42 ergs s" 1 



0.768±0.086 



(1) 

The intrinsic scatter about this relationship is w 0.179 
dex, implying a potential uncertainty in ljjy inferred 
from this relationship of the same magnitude. There was 
no trend in the residuals with either z or Lh/3, implying 
that Equation fl} should give unbiased estimates of Ijjv 
for the z < 0.4 sources. Values of Ljjv were estimated 
using Equation (P) for both the K07 and S05 z < 0.4 
sources, a total of 42 sources. 

The distributions of Ljjy and Lx as a function of red- 
shift are shown in Fi gure [H We calculate t he ratio of 
optical to X-ray flux ( Tananbaum et aLl ll979) as 

(2) 



a ox = -0.384 log(/ 2kc v//250o), 



where /2koV and /2500 are the rest-frame flux densities 
at 2 keV and 2500A, respectively. If the flux density 
from 2500A to 2 keV is a simple power law, then a ax is 
the spectral slope of this continuum, and thus a ox may 
be thought of as a crude estimate of the shape of the 
ionizing continuum. The parameter a ox is an impor- 
tant parameter for model comparison, as it summarizes 
the amount of energy emitted in the X-ray region (most 
likely a Comptonizcd component), compared with that 
emitted in the optical- UV (accretion disk component). 
The distribution of function of Ljjy and z are 

also shown in Figure [T] 

The error on a ox is the result of measurement errors on 
the UV and X-ray flux, as well as error caused by quasar 
variability over the different epochs for the UV and X- 
ray observations. In general, the measurement errors on 
luv are negligible compared to the error on lx- S05 esti- 
mates a error on the X-ray flux of ~ 0.23 dex, including 
both the contributions from measurement error and vari- 
ability. Typical long-term X-ray variability for Seyfert 
Is is 20%-40% with no obvious trend with luminosity 
dGrupe et al.ll200U lUttlev et al.ll200l iMarkowitz et all 
2003). The measurement errors in lx for the K07 sam- 
ple are typically ~ 0.07 dex. Assuming X-ray variability 
amplitudes of 30%, this implies typical uncertainties in 
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the X-ray luminosity of ~ 0.15 dex. Therefore, we esti- 
mate the uncertainty on a x to be ~ 0.06 for the K07 
sources and ~ 0.09 for the S05 sources. 

3.2. Line Profile Extraction and Fitting 

We extracted the H/3, Mg II, and C IV emission 
lines in order to u se their widths in our black hole 
mass estimates (e.g.,IVe stcr£aard 2002; M cLure fe Jarvisl 
I2002t IVestergaard fe Petersonll2006l ) . These line were ex- 
tracted by first subtracting the continuum and Fe emis- 
sion, interpolating over any narrow absorption features, 
and modelling all lines within the extraction region as a 
sum of Gaussian functions. Any nearby lines were then 
subtracted, leaving only the broad emission line profile. 
In all cases the line profile extraction was done interac- 
tively and every line fit was inspected visually. 

For H/3, we extracted the region within ±2 x 
10 4 km s _1 of 4861 A, where we use the standard con- 
vention that negative velocities are blueward of a given 
wavelength. The H/3 profile was modeled as a sum of 
2-3 Gaussian functions. The [O III] A4959A and [O III] 
A5007A lines were modeled as a sum of 1-2 Gaussian 
functions, depending on the signal-to-noise of the lines. 
A sum of two Gaussian functions was used for the higher 
S/N lines because the [O III] line profiles are not ex- 
actly a Gaussian function; the individual Gaussian com- 
ponents are not considered to be physically distinct com- 
ponents. The widths of the narrow Gaussian functions 
for H/3 and [O III] lines were fixed to be equal to ea- 
chother. The [O III] lines and the narrow component 
of the H/3 line were then subtracted, leaving the broad 
component of H/3. 

For Mg II, we extracted the region within ±2 x 
10 4 km s _1 of 2800A. There are no nearby non-iron emis- 
sion lines that Mg II is blended with, so the extraction 
is trivial after removing the Fe and continuum emission. 

For C IV, we extracted the region within —2 x 
10 4 km s~ x and 3 x 10 4 km s" 1 of 1549A. The C IV 
line was modeled as a sum of 2-3 Gaussian functions, 
and He II A1640 and O III] A1665 were modeled as a 
sum of 1-2 Gaussian functions each. After obtaining es- 
timates of the He II and O III] profiles, we subtracted 
these components. We did not model the N IV] A1486 
emission line as this line is typically weak and lost in the 
C IV wings. 

In order to estimate Mbh, it is necessary to measure 
the FWHM of the emission lines. After extracting the 
line profiles, we estimate the FWHM for the H/3, Mg 
II, and C IV emission lines by fitting them to a sum of 
1-5 Gaussian functions, enabling us to obtain a smooth 
representation of each line. In contrast to our profile ex- 
traction technique, we choose the number of Gaussian 
functions t o minimize the Bayesian Information Crite- 
rion (BIC, lSchwartall979f ). The BIC is a common cri- 
terion to use for selecting the number of parameters i n 
a model (e.g., see lHastie. Tibshirani. fe Friedmanll200lD ; 
the model that minimizes the BIC is approximately the 
model that is most supported by the data. For Gaus- 
sian errors, as assumed in this work, the BIC is simply a 
modification to the standard x 2 statistic: 

BIC = x 2 + 3Klnn, (3) 

where K is the number of Gaussian functions used, 3K 
is the number of free parameters, and n is the number 



of data points used in the fit. Using the BIC to 'fit' 
the number of Gaussian functions thus allows us more 
flexibility in obtaining a smooth representation of the 
line profile, as we are not choosing the number of Gaus- 
sian functions arbitrarily. Once a smooth representation 
is obtained, we automatically measure the FWHM di- 
rectly from the best fit line profile. 

The standard errors on FWHM are estimated using a 
bootstrap method. We simulated 100 'observed' emission 
lines by adding random Gaussian noise to the best fit line 
profile with standard deviation equal to the noise level 
of the spectrum, including the propagated errors from 
the continuum and iron emission fitting. We then fit 
each of the simulated emission lines, keeping the number 
of Gaussian functions fixed at the number found from 
fitting the original profile, and measured the FWHM for 
each simulated line. The standard error on FWHM was 
then estimated as the standard deviation of the FWHM 
values measured from the simulated line profiles. 

4. ESTIMATING M BH 

Recently, reverberation mappin g studies of broad line 
AGN (e.g., iPeterson et al.ll200l established a correla- 
tion between the broad line region (BLR) size, R, and 
the continuum luminosity (the R-L relationship, e.g., 
iKaspi et al] l2005t iBentz et all I2006D . This has made 
it possible to estimate black hole virial mass Mbh = 
fv 2 R/G for individual sources, where the BLR velocity 
v is estimated from the width of an emission line (e.g., 
i Wandel et aLlll999t IVestergaard! I20 02t iMcLure fe Jarvisl 
l2002t IVestergaard fe Petersonir2006l ) . We choose the pro- 
portionality constant to give broad line mass estimates 
consi s tent with the Mbh-o ~ relat i onship ([G cbhard t~et al.l 
I2000t iMerritt fe Ferraresd [200ll: iTremaine et al.l 12002). 
/ = 1.4 ± 0.45 (|Onken et all l2004f) . An estimate of 
the Eddington luminosity can be computed as LEdd — 
1.3 x 10 3S M B h/Mq erg s" 1 . 

4.1. Black Hole Mass Estimates from H[3, Mg II, and C 

IV 

In this work we estimate Mbh from the H/3, Mg II, 
and C IV emission line s . W e use the relationship of 
IVestergaard fe Petersonl (|2006h to estimate Mbh from 
the H/3 and C IV emission lines, and the datails of the 
Mg II calibration will be dis cussed in a forthcoming pa- 
per IVestergaard et all (|2008l in preparation) . The cali- 
bration for the Mg II mass estimates was calculated to 
ensure that they are consistent with the mass estimates 
based on H/3 and C IV. We have 49 sources with both H/3 
and Mg II mass estimates, and 73 sources with both C 
IV and Mg II mass estimates. Both samples show consis- 
tent mass estimates between the different emission lines, 
within the intrinsic uncertainty in the broad line mass 
estimates (~ 0.4 dex). 

We will denote the broad line mass estimates as Mb l , 
and rriBL = log Mbl/Mq. It is important to distinguish 
between M B l and M B h, as M BL oc L^FWHM 2 , 7 ~ 
0.5, is an estimate of Mbh derived from reverberation 
mapping, and thus in general Mbl 7^ Mbh- The 
statistical uncertainty needs to be taken into account 
when analyzing correlations involving derived quantities 
like Mbl L^FWHM 2 , as they can bias the results 
(|Kellv fe BechtoldlM7tlK^l200l . 
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The uncertainty in / increases the formal statistical 
uncertainty in the broad line estimates of Mbh to ~ 0.46 
dex. Our adopted formal uncertainty of ~ 0.46 dex is 
merely statistical, and additional systematic unc ertain- 
ties in reverberation m apping may contribute (jKrolikl 
[2001 lCollinlt~aT][200l . Because our adopted uncer- 
tainty of ~ 0.46 describes the scatter in broad line mass 
estimates about the reverberation mapping estimates, as 
calibrated via the Mbh~o~ relationship, the regression re- 
sults found in this work should be understood as results 
that could have been obtained if we had reverberation- 
based Mbh for the sources in this work. However, in- 
stead of reverberation-based mass estimates, we have 
broad line mass estimates with 'measurement error' equal 
to ~ 0.46 dex with respect to the reverberation-based 
mass estimates, thus increasing the uncertainty from the 
regression analysis. 

For most sources, measurement errors on FWHM and 
L\ did not significantly contribute to the uncertainty on 
Mbh- If there were two emission lines in the same spec- 
trum we averaged the two mass estimates, where the av- 
erage was weighted by the uncertainties in the two esti- 
mates. The distribution of Mbh as a function of z for 
our sample is shown in Figure [5] The broad line masss 
estimates for the sources in our sample are reported in 
Table HI 



tive approach and merely compare a ox and T x with 
Ljjv I 'LEdd and Lx/LEdd- The estimated values of 
Luv I LEdd and Lx/LEdd for the sources in our sam- 
ple are reported in Table [TJ We can write Ljjv I LEdd oc 
fuvLboi/ LEdd and Lx/ LEdd oc fxLboi/ Lsdd, where fuv 
and fx are the inverses of the bolometric corrections for 
Luv and Lx, respectively. The quantities fuv and fx 
are proportional to the fraction of the bolometric lumi- 
nosity emitted at 2500A and 2 keV. Then, correlations 
between either a ox or Tx and Luv / LEdd will result if 
a ox or is correlated with fuv, the Eddington ratio, 
or both, and likewise for fx- 

While we do not use an estimate of the Eddington ratio 
in our analysis, it is helpful to estimate the distribution of 
Eddington ratios probed by our sa mple. We assume the 
bolometric correction described in iHopkins et alj (|2007f ) 
for the z < 1.5 sources, an d constant bolometric de- 
scribed in lVestergaardl (|2004fl of L bol = 4.62A J L A (1350A) 
at z > 1.5. In Figure [3 we also show the distribution of 
estimated Eddington ratios as a function of z. Because 
the distribution of estimated Lboi/ L Edd is the true dis- 
tribution of Lboi I L Edd broadened by the distribution of 
errors in the estimates, our sample likely probes a smaller 
range in Eddington ratio than that inferred from Figure 
[2j Therefore, at most our sample probes RQQs with 
Eddington ratios 0.03 < Lboi/ LEdd 2. 



4.2. Eddington Ratio Estimates 

A constant bolometric correction has been used in most 
previous studies involving the AGN Eddington ratio . 
However, recent work by IVasudevan &: Fabian! (|2007l ) 
has suggested that bolometric corrections show a large 
spread with no obvious dependence on luminosity. Fur- 
thermore, these authors found evidence that the bolo- 
metric correction depends on the Eddington ratio. This 
implies that the error in the bolometric correction is cor- 
related with Eddington ratio, which therefore implies 
that the error in the estimated Eddington ratio is cor- 
related with the actual Eddington ratio. An Eddington 
ratio-dependent error in the bolometric correction my 
cause problems when using the estimated Eddington ra- 
tios to infer correlations. 

Further difficulties with a constant bolometric correc- 
tion are illustrated with Figure [3] In [3] we plot a ox as 

a function of Luv / LEdd and Lx/LEdd, where we esti- 
mate the Eddington luminosity from the broad line mass 
estimates as L E dd = 1-3 x 10 38 M B l/M q erg s _1 . As 
with Mbl, we use the notation LEdd to emphasize that 
L Edd is an estimate of the true LEdd based on the broad 
line mass estimates, and therefore Lsdd oc L 1 FWHM 2 . 
Constant bolometric corrections are often applied to ei- 
ther the optical/UV or X-ray luminosity. If a constant 
bolometric correction was valid for both Luv and Lx, 
then we would expect that Luv / LEdd oc Lx/LEdd oc 
Lboi I LEdd- However, while a correlation between a ox 
and both Luv / LEdd and Lx/LEdd is apparent, they 
are of opposite sign. Because the correlations are of 
opposite sign, it cannot be true that both Luv I LEdd 
and Lx/LEdd are proportional to the Eddington ratio, 

Lbol I 'LEdd- 

Because of the current significant uncertainty regard- 
ing RQQ bolometric corrections, we take the conserva- 



5. DEPENDENCE OF a Q x ON M BH ,L UV /L EDD , AND 
Lx/Ledd 

We used our sample of 318 sources with estimates of 
Mbh to investigate the dependence of a ox at a given 
black hole mass, Luv / LEdd, and Lx/LEdd- We use lin- 
ear regression analysis in order to understand how a ox 
varies wit h respect to these parameters. We use the 
method of iKellvi (|2007f ) to est i mate the regression pa- 
rameters. The method of lKellvi (|2007f) accounts for mea- 
surement er rors, non-det ections, and intrinsic scatter. 
In addition, IKellvi (|2007f ) adopts a Bayesian approach, 
computing the posterior probability distribution of the 
parameters, given the observed data. Thus the uncer- 
tainties on the regression coefficients have a straight- 
forward interpretion, and do not rely on large-sample 
approximations. Many other methods, such as tradi- 
tional maximum-likelihood, assume that the errors in 
the regression parameters follow a Gaussian distribu- 
tion, which is valid as the sample size approaches infinity. 
However, this assumption is not necessarily valid for our 
finite sample size, especially in the presence of censor- 
ing (i.e., presence of upper/lower l imits) and si gnificant 
measurement error. The method of lKellvi (|2007f ) directly 
estimates the probability distribution of the regression 
parameters, and is therefore preferred. 

We assess the simple 2-dimensional correlations be- 
tween a ox and M B h, L V v / L E dd, and L x /L E dd, and 
compare with the a ox -Luv correlation. The results from 
the regressions are 

a ox = -3.91+|« + (0.12t°;°2) log L uv , o am = 0.14+°;°?, 

p = 0.57+°;°° (4) 

a ox = 0.05±g;l + (0.17+°;°!) log M BH , cr Qox = 0.14+°;°J, 

p = 0.53+°;^ (5) 

a ox = 2.90tg;S + (0.99±g;|S) log L uv /L E dd, <r aox = O.Q5+°;°|, 
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o.95±g:2S 

h0.25 
-0.34 



- (0.57+°;™) log L x /L Edd , a aox = 0.03 



(6M2500A) _ +1 . 02 ( vL v (2 keV) 
-0.04 r (2keV)- U - 85 -°- 52 l, L Edd 



0.03' 



-1.48±0.14 



(11) 



p = -o.98ig:g 

where <7 Qox is the intrinsic dispersion in a ox at a given 
Luv,M BH ,L uv /L Edd , or L x /L Eddl p is the linear cor- 
relation coefficient for a ox and the respective indepen- 
dent variables, and the errors are quoted at the 95% (2a) 
level. All four relationships are significant, with RQQs 
becoming more X-ray quiet as Ljjv, M E h, or Ljjy I L Edd 
increases, and more X-ray loud as Lx/L Edd increases. 
Because we have attempted to correct for the intrinsic 
statistical scatter in the broad line mass estimates, Equa- 
tions ([4])-([7]) refer to the intrinsic relationships involv- 
ing Mbh, barring any systematic errors in reverberation 
mapping, and are not simply correlations between a ox 
and the broad line mass estimates. The estimated dis- 
tributions of a ox as a function of Ljjv, M E h, Ljjv / L Edd , 
and L*x / L Edd are shown in Figure [4j along with the re- 
gression results. 

The intrinsic dispersion in a ox quantifies the magni- 
tude of scatter in a ox at a given Luv, M E h, Ljjv / L Edd , 
or Lx/L Edd . Because we have attempted to account for 
contribution to the scatter in a ox resulting from measure- 
ment error and variability, <7 Qox represents the dispersion 
in the real physical scatter in a ox over the population of 
RQQs. This 'residual' scatter represents the amount of 
variation in a ox that is not 'explained' by variations in 
L uv ,M B HiL uv /L Edd , or L x /L Edd , respectively. This 
intrinsic scatter in a ox may be due to variations in accre- 
tion rate, viscosity, column density, and other quantities 
not included in our regression. 

In §[2] we estimate that there are 25 ± 5 BAL quasars 
in our sample at z < 1.5. Because these objects have 
the appearance of being X-ray weak, and because red- 
shift is artificially correlated with luminosity and Mbh 
in a flux limited sample, we expect that the presence of 
unidentified BALs at z < 1.5 will produce an excess of 
X-ray weak objects at low Ljjv and Mbh, thus flatten- 
ing the inferred slopes. Inspection of the plot of a ox and 
z in Figure [1] suggests an excess of X-ray weak objects 
at z < 1.5 and a x ^ 1-8, implying these objects are 
BAL quasars. We removed these 10 objects and refit the 
regressions. Omission of these objects resulted in a steep- 
ening of the slopes for the Luv and M E h regression, and 
a flattening of the slope for the Luv / L Edd regression. In 
addition, the intrinsic dispersion in a ox decreased for the 
Luv and M E h regressions, while it remained the same 
for the Ljjv / L Edd regression. These changes were small 
(~ 10%) and have no effect on our conclusions. There 
was no difference in the results for the Lx / L Edd regres- 
sion. 

Once can use Equation ([2]) to express the regression 
results (Eq. [4] [7]) in the alternate form 



M2500A) +0 . 08 4 ^(2500A) ^ 

L„(2 keV) -^'-o-or xw { IC^ergs" 1 J 

/ , , \ 0.43±0.06 
I +0.65 „ n 3 / M BH 



0.31±0.03 



M2500A)_ , 
L„(2keV) - J - 8i -°-63 x Vi O 9M 

L.^OOA) 

L v (2 keV) ~ 2 - jS V 



(8) 
(9) 



2.57±0.45 



7(>'LMrm\)Y ( . |(|) 



(jfjheve the intrinsic dispersion in log Ljjv /Lx at 
a given L uv , M B h, L uv / L Edd , and L x /L Edd is ~ 
0.356,0.375,0.133, and 0.089 dex, respectively. In con- 
trast to Equations ([4])-([7|), we quote the 68% (la) un- 
certainties on the constants of proportionality, and the 
posterior standard deviations on the exponents. Equa- 
tions (f8"|)- (fTTj) may be more physically interpretable and 
allow easier comparison with models. 

6. NONMONOTONIC DEPENDENCE OF T X ,L UV /L EDD , 
AND L x /L B ol 

Recent work has suggested a correlation between 
quasar X-ray spectral slope, ax — — oc 
v~ ax , and quasar Eddington ratio as inferred from 
broad line mass estimates based on the H/3 emission 
line (e.g., iPorauet et al.l l2004t iPiconcelli et all 120051 : 
IShemmer et all l2006[ )~The iKellv et all (|200'7f ) sample 
contains measurements of Tx for 157 sources, and we 
were able to estimate black hole masses for 153 of them. 
In this section we use these 153 RQQs to investigate the 
dependence of T x on M B h,L uv /L Edd , and L x /L Edd . 

6.1. Regression Analysis 

The distributions of Tx as a function of estimated 
black hole mass, Ljjv / L Edd , and Lx/L Edd are shown in 
Figure[H]for the entire sample, and in Figure [7] seperately 
for each emission line. While there does not appear to be 
a monotonic trend between Tx and M E h, Ljjv / L Edd , or 
Lx/L Edd when using the entire sample, there is evidence 
for a trend between Tx and these quantities when using 
the H/3 line, and an opposite trend between Tx and these 
quantities when using the C IV line. 

We performed a linear regression of Tx on 
log M B h , log L uv / L Edd , and log L x / L Edd seperately for 
each e missio n line. As before, we used the method of 
iKellvl (|2007f ) when performing the regression in order 
to correct for the intrinsic statistical uncertainty in the 
broad line estimates of Mbh- The results for M E h are 

T x = 5.69±£.29 " (°- 44 ±o.9i) log Mbh, ° = 0-43±g;i?, 

p = -0.4518:2, (H/3) (12) 
T x = -19.01&8- (l.92l£S) logM BH , <7 = 0.27i{|;& 
p = -0.621^, (Mgll) (13) 
(Q.52iS;£) lo E M BH , a = 0.22^.11, 



^ = -2.79t? 22 8 



0.53: 



0.58' 



(CIV) 



(14) 



J Edd 



the results for Ljjv I 'L Edd are 

T x =ZM + _ 2 1 f 5 + (l.23tS;$|) lo S L UV /L E dd, <r = 0.26lg;g, 

p = 0.87tg;& (H/3) (15) 

T x = 5. 13tll 8 4 + (2.UtHl) log L uv /L Edd , ° = 0-^t\]f 2 , 

p = 0.541?^, (Mgll) (16) 

r x = 0.85l°;^ - (0.95±?;||) logL uv /L Edd , a = 0.17±° :{l 

p = -0.8118;™, (CIV) (17) 

and the results for Lx/L Edd are 

T x = 4.24l 2 ;g + (0.85lJ;°J) log L x /L Edd , a - 0.33+°i?, 
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P = 0.761°;^, (H/3) (18) 
T x = 7.26™ + (1.97^) log L x /L Eddl a = 0-23±°:i|, 

P = 0.77l°; 2 ^, (Mgll) (19) 
r x = -0.54±£gj _ (o.96±| : 5|) log L x /L Edd , a = 0.21 i^. 

P=-0.65+ J°(CIV). (20) 



In these equations we have quoted the errors at 95% (2a) 
confidence. The probability distributions of the slope 
and intrinsic dispersion are shown in Figure[5] The larger 
uncertainty in the results for the Mg II sample is likely 
caused by the more narrow range in Luv, L x , and Mbh 
probed. 

There are formally no significant linear correlations 
for the T x -Mbh relationship. However, there is a sta- 
tistically significant difference between the H/3 and C 
IV slopes, with w 99.3% of the posterior probability at 
Prn > 0rn where (3 m denotes the T x -\ogMsH regres- 
sion slope. The probability distribution for the difference 
in slopes from the Mbh regression is shown in Figure [5J 
The significant difference in the slope for the H/3 and C 
IV sample implies a nonlinear relationship between T x 
and ttibHi in spite of the fact that the H/3 and C IV corre- 
lations themselves are not 'statistically significant'. Re- 
sults similar to the T x -\ogMBH regressions were found 
for the T x -Ljjy / Lbm and T x -L x / Ljjy regressions, but 
with opposite sign and higher statistical significance. 

We performed monte carlo simulations as a consistency 
check on our inferred non-monotonicity of the T x rela- 
tionships. While we have attempted to account for the 
significant statistical uncertainty on the broad line mass 
estimates, we employ these monte carlo simulations to 
ensure that the observed non-monotonic behavior is not 
a spurious result caused by the uncertainty on Mbh- Wc 
performed 10 5 simulations under two null hypotheses: 
(1) that T x is independent of L\jy / 'L Edd , and (2) that 
T x depends linearly on log Luy/L Ed d- F° r both cases 
we simulated black hole mass estimates derived from H/3 
and C IV seperately. We first simulated 'true' values 
of Mbh for each emission line from a normal distribu- 
tion with means equal to the observed mean of the two 
respective subsamples, and variances equal to the differ- 
ence between the observed variance of the subsamples 
and the average intrinsic variance in the broad line mass 
estimates. To simulate the uncertainty in the mass esti- 
mates, we added random Gaussian errors to these 'true' 
values of Mbh with standard deviation equal to the un- 
certainty in the mass estimates, ~ 0.4 dex. For the case 
where T x was assumed to be independent of Ljjy / L E dd, 
we simulated values of T x from a normal distribution 
with mean equal to the sample mean of T x and variance 
equal to the difference between the observed variance in 
T x and the average of the variance in the measurement 
errors. For the case where T x was assumed to depend 
linearly on log Lyy / 'L E dd, we simulated values of Tx ac- 
cording to our best fit relationship to the H/3 subsample, 
given by Equation (fT5| . Finally, for both cases we added 
random Gaussian errors to the simulated values of T x by 
randomly reshuffling the dispersions in the measurement 
errors in T x . 

For each of the 10 5 simulated samples, we selected 
those samples that displayed a non-monotonic trend, i.e., 
those sample where the slope for the H/3-based regression 



had a different sign from the slope of the C IV-based re- 
gression. Under the hypothesis that T x is independent 
of L\jy /L E dd, only 3 of the I0 5 simulated samples had 
both a non-monotonic trend and an absolute value of 
the difference in slopes between the H/3 and C IV regres- 
sion that were larger than that observed for our actual 
sample. Under the hypothesis that T x depends linearly 
on log Ljjv / L E dd, none of the 10 5 simulated samples ex- 
hibited a non-monotonic trend. Therefore, the observed 
non-monotonic trend in T x with Eddington ratio is not 
a spurious result caused by the statistical uncertainty in 
the broad line mass estimates, in agreement with our 
Bayesian regression results. 

In order to investigate whether the non-monotonicity 
in the dependence of T x on Eddington ratio depends on 
Mbh, we performed a linear regression of T x simulta- 
neously on log Luy I 'L E dd and log Mbh- This also al- 
lows us to quantify whether the Eddington ratio is the 
driver behind the T x -M E h relationship. In particular, 
the T x -Mbh relationship is weak compared to the Ed- 
dington ratio relationships, and therefore it is reasonable 
to conclude that Eddington ratio is the primary driver in 
these relatio nships . We applied the multiple regression 
technique of lKellvl (|2007[ ) seperately to both the H/3 and 
C IV subsamples. The results are: 

rjc = 2.58±t|S + (0.18t°j°) log Mbh + 

(l.32tS$) log WW, o = 0.28t° 'j2, (H®1) 
T x = -4.26+^ + (OMtlH) log M BH - 

(0M±°f 6 )logL uv /L Edd , a = 0.l7t° ii (<TO 

Here, we have quoted the errors at 95% significance. 
There is no statistically significant evidence that T x de- 
pends on Mbh at a given Ljjy / L E dd, and therefore we 
conclude that the primary driver in the T x relationships 
is Eddington ratio. However, this does not rule out the 
possibility that the non-monotonic trends with Edding- 
ton ratio are the result of a discontinuous change in the 
slopes at a 'critical' Mbh, as discussed in the next two 
sections. 

6.2. Is the Sign Change in the Correlations Caused by 
the Different Emission Lines Used to Estimate 
M BH ? 

The opposite correlations for H/3 and C IV are intrigu- 
ing but may represent problems with the broad line mass 
estimates. In particular, it is possible that the error in 
the broad line mass estimates is correlated with T x , but 
in opposite ways for H/3 and C IV. The most likely source 
of such a spurious correlation would be a correlation be- 
tween T x and the scatter about the R-L relationship 
for H/3 and C IV, respectively. For example, if one were 
to systematically overestimate R with increasing T x for 
the C IV emitting region, then one would infer a larger 
Mbh from C IV, and thus one would infer a spurious 
correlation between Mbh and T x . However, the H/3 line 
is only available at z < 0.8 and the C IV line is only 
available at z > 1.6, and thus the change in sign for the 
T X -M E H correlations could be due to different spectral 
components shifting into or out of the observable X-ray 
spectral region (0.3-7.0 keV). In Figure IT0| we show the 
distribution of Mbh as a function of z for the H/3 and C 
IV samples. As is clear from Figure [101 the C IV line is 
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probing sources with Mbh ^ 3 x 10 8 Mq, while the H/3 
line is probing sources with Mbh 3 x lO 8 Ai0. There- 
fore, the change in sign for the Tx correlations could also 
be due to something more physically interesting, such as 
a change in the structure of the corona that occurs at 
some critical black hole mass or accretion rate. 

We can test if the change in sign for the Tx correla- 
tions is the result of problems with the mass estimates for 
either H/3 or C IV, or if it is the result of the differences 
in z and Mbh probed by the two lines. While the I\x— 
Mbh relationship is weak compared to the Eddington 
ratio dependencies, we can use the Tx~Mbh relation- 
ship to test whether the non-monotonic dependency of 
on Eddington ratio is a spurious result caused by 
systematic difference in the broad line mass estimates. 
This is because the Eddington ratios are inferred from 
the broad line mass estimates, and any systematic differ- 
ences between the Eddington ratio inferred from H/3 as 
compared to C IV should also manifest themselves in the 
weaker r x~Mbh relationship. This is true irregardless 
of whether the difference in slopes between the H/3 and 
C IV Tx-Mbh relationships is 'statistically significant' 
or not. 

We compiled five sources from the literature at z > 1.3 
with H/3-based mass estimates of Mbh > 5 x 10 9 M q . 
In addition, we compiled five more sources from the lit- 
erature at z < 0.2 and C IV-based mass estimates of 
Mbh < 10 8 M Q . This 'test sample' of 10 sources is listed 
in Table [21 Then, we test whether the H/3 test sources 
are better described by the H/3 regression or by the C IV 
regression, and likewise for the C IV test sources. If the 
change in sign for the Tx correlations is due to problems 
with the broad line mass estimates, then we would expect 
the H/3-based mass estimates to be better described by 
the H/3 regression. However, if the change in sign is due 
to the difference in redshift and Mbh probed by the two 
regressions, then we would expect the H/3 test sources to 
be better described by the C IV regression, as the H/3 
test sources are at high-z and have high-Ms//. A similar 
argument applies to the C IV test sources, since they are 
at low-z and have low black hole masses. 

Figure [Til compares Tx and Mbh for H/3 and C IV for 
both the sources in our main sample and the test sources, 
as well as the best fit regression lines for the H/3 and 
C IV samples, respectively. The high-z, high-Ms^ H/3 
test sources appear to be better described by the high-z, 
high- Mbh C IV-based regression, and likewise the low-z, 
\ow-Mbh C IV test sources appear to be better described 
by the low-z, Iow-Mbh H/3-based regression. 

We can quantify this result by calculating the proba- 
bility that the H/3 test sources 'belong' to the H/3-based 
regression, as compared to the probability that the H/3 
test sources 'belong' to the C IV-based regression. As- 
suming that the test sources are as equally likely to be- 
long to either regression a priori, this ratio of proba- 
bilities is simply the ratio of the likelihood functions of 
the test sources for each regression relationship, where 
the likelihoo d functions are given by Equation (24) in 
Kellvl (120071); t his ratio is called the 'Bayes Factor' (e.g., 
CongdonH 2006) . In order to incorporate our uncertainty 
in the regression parameters, we use the value of the 
likelihood function averaged over the probability distri- 
bution of the regression parameters. We find that the H/3 
test sources are « 250 times more likely to 'belong' to 



the C IV-based regression, and that the C IV test sources 
are « 140 times more likely to 'belong' to the H/3-based 
regression. Because the test sources are independent, it 
follows that the test sources are > 10 4 times more likely 
to be described by the regression fit using the opposite 
emission line sample. This is further evidence that Tx 
and Mbh are not statistically independent; if Tx and 
Mbh were independent, then the test sources would not 
show a strong preference for either regression. Based on 
this analysis, we conclude that the change in sign of the 
Tx correlations is not due to problems associated with 
the use of the H/3 and C IV emission lines, but rather 
due to the different range of z and Mbh probed by the 
two subsamples. 

6.3. Is the Sign Change in the Correlations Caused by 
the Different Redshift Ranges Probed? 

While it appears that Tx has a nonmonotonic depen- 
dence on Mbh, Luv /L<Edd, and Lx/LEdd, it is unclear 
as to whether the sign change in the correlations is depen- 
dent on z or Mbh- However, we can test this in the same 
manner as was used to test if the sign change is due to 
problems with the H/3- or C IV-based mass estimates. In 
this case, we need a sample of high-z, \ow-Mbh sources 
in order to break the degeneracy between Mbh and z. If 
the sign change in the correlation is redshift dependent, 
then we would expect the test sources to be better de- 
scribed by the high-z, high-MsH C IV-based regression; 
but, if the sign change is black hole mass dependent, then 
we would expect the test sources to be better described 
by the low-z, \ow-Mb h H/3-based regression. 

Our test sample consists of nine z > 1,Mbl < 3 x 
10 8 M (7) quasars from the Cosmic Evolution Survey (COS- 
MO S. IScovilld 12007ft. w ith optical spectra from Magel- 
lan (iTrump et al.l [2007) and X-ray spectra from XMM- 
Newton (Mainicri 2007). Black hole masses for these ob- 
jects ([Trump et al.ll2008l . in preparation) were estimated 
from the Mg II and C IV emission lines in the same 
manner as above. The objects are summarized in Table 
[3j and their location in the Mbh~z plane are shown in 
Figure [T2l The COSMOS sources break the degeneracy 
between Mbh and z present in our SDSS sample, and 
are therefore adequate to test for a redshift dependence 
in the slope of the Tx~Mbh relationship. In Figure [13] 
we compare the COSMOS test sources with the H/3- and 
C IV-based regressions. As can be seen, the high-z, low- 
Mbh COSMOS sources are better described by the low- 
z, Iow-Mbh regression. We can quantify this in the same 
manner as described in § 16.21 by averaging the likelihood 
function of the test sources over the posterior probabil- 
ity distribution. We find that the COSMOS test sources 
are > 10 5 times more likely to be better described by 
the TL[3-T x regression, and thus the sign change in the 
Tx correlations is not due to the difference in redshifts 
probed by the H/3 and C IV samples. 

7. DISCUSSION 

Previous work has found e yidence for a corre lation 
between a a ^ and both Mbh (|Brunner et al.lll997ft and 
Lboi/LEdd ([Wang et al.ll2004ft . and for a correlation be- 
tween T x and both M B h (e.g., iPorquet et all 120041: 
Piconcelli et alj l2005ft and Lboi/LEdd (e.g., )Lu k, Yul 
1999: iGierlihski fc Pond 12003: [Shcmmer et al. 2006ft . in 
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aggreement with the results found in this work. How- 
ever, our study differs from previous work in that we 
study a large sample of RQQs (318 sources with a ox , 
153 with Tx) over a broad range in black hole mass 
(10 6 < Mbh/Mq < 10 10 ) and redshift < z < 4.8; 
to date, this is the largest study of the dependence of 
the X-ray properties of RQQs on Mbh, Ljjv / LEdd, and 
Lx/LEdd- In addition, this the first study of its kind to 
correct for the intrinsic statistical uncertainty in broad 
line mass estimates when quantifying the intrinsic trends 
between the X-ray emission and Mbh, Ljjv / LEdd, and 

Lx/LEdd- 

Currently, there are two main types of geome- 
tries being considered for the comptonizing corona. 
The first of these is that of a 'slab '-type geome- 
try, possibly patchy, that sandwiches the disk (e.g. , 
Bisnovatvi-Kogan fc Blinnikovlll977t iGaleev et alj|1979t 
Navakshinl 120001 : ISobolewska et al.ll2004bf ). and the sec- 
ond is that of a hot spherical inner advection dom- 
inate d flow (e.g., IShapiro et al.l Il976t [Zdziarsk i et al.l 
1999); hybrids b etween the two geometries have also been 
consid ered (e.g.. iPoutanen et al.lll997t ISo bolewsk a et all 
[20n4ah . There is a growing body of evidence that the ad- 
vection dominated hot inner flow does not exist in objects 
with Eddington ratios L bo i/ LEdd ^ 0.01, as inferred from 
the existence of a r e lativistica l ly broadened iron line (e.g . , 
iMineo et all 120001 : iLee et al.l 12001 iFabian et alll2002h . 
relativistically broa dened reflection of ionized material 
(jJaniuk et al.ll200lh. and by analogy with galact ic black 
holes (e.g., lEsin et all Il997t iNowak et alll2002l ). The 
range in Eddington ratios probed by our study is at most 
0.03 < L bo i/L E dd < 2, with a mean of L bo i/L E dd ~ 0.25. 
Therefore, the RQQs in our study are likely to have disks 
that extend approximately down to the last marginally 
stable orbit, and thus should only have the 'slab' type 
geometries. 



7.1. Dependence of a ox on Mbh 

In this work we have found that RQQs become more 
X-ray quiet as Mbh increases, and confirmed the well- 
established relationship between a ox and Ljjv- Because 
Ljjv increases with Mbh and the accretion rate relative 
to Eddington, to, the well-known a ox -Ljjy correlation is 
likely driven by the C( ox -Mbh and a ox -rh correlations. 
A correlation between a ox and Mbh is expected even if 
the fraction of the bolometric luminosity emitted by the 
disk is independent of Mbh, as the effective temperature 
of the disk depends on Mbh- As Mbh increases, the 
effective temperature of the disk decreases, thus shifting 
the peak of the disk emission toward longer wavelengths. 
Because the flux density at 2500A lies redward of the 
peak in the disk SED over most of the range Mbh probed 
by our study, this shift in the disk SED toward longer 
wavelengths produces an increase in Ljjv relative to Lx ■ 

We can use the standard thin disk solution to assess 
the evidence that the fraction of energy emitted by the 
corona depends on Mbh- We assume a simple model 
where the spectrum for the disk emission is that expected 
for an extended thin accretion disk, and the spectrum for 
the corona emission is a simple power-law with exponen- 
tial cutoffs at t he low and high energy end. According to 
Wandel (2000), the spectrum from a radially extended 



thin accretion disk can be approximated as 

finA D (— ) V"/"-, (23) 

where Ad is the normalization and v co is the cut-off fre- 
quency. In this work, we choose the normalization to 
ensure that Equation (|23p integrates to unity, and there- 
fore /jP gi ves the shape o f the disk emission. For a Kerr 
black hole, Malkan (1991) finds that the cut-off frequency 
is related to Mbh as 

hv co = (6eV)rh 1 / 4 (M S K/lO 8 M )- 1 / 4 . (24) 

We assume that the X-ray emission from the corona can 
be described by a simple power law with an exponential 
cutoff at the high and low end: 

fC = A c v- {Tx -^e- v/Uh ^ h e- Vl °™l v . (25) 

Here, Ac is the corona spectrum normalization, v^igh is 
the high energy cutoff, and i>i ovl is the low energy cutoff. 
We choose the low energy cutoff to be ui ow = 20 eV, and 
we c hoose the high en ergy cutoff to be Vhigh — 200 keV 
(e.g.- lGilli et al.|[2007t ). As with Equation ([23]) . we choose 
the normalization in Equation (|25j) to be equal to unity. 

Denoting fp to be the fraction of bolometric luminos- 
ity emitted by the disk, our model RQQ spectrum is then 

Lu « L bol [f D f? + (1 - f D )tf] . (26) 

We computed Equation ([26]) assuming a value of rh = 0.2 
and fr> = 0.85. We chose the value of the to = 0.2 be- 
cause it is representative of the RQQs in our sample, and 
we chose the value fo = 0.85 because it gives values of 
a ox typical of the RQQs in our sample. We vary Mbh 
but keep fjj and rh constant because we are interested 
in investigating whether there is evidence that assum- 
ing independence between Mbh and both fjj and rh is 
inconsistent with our a ox results. 

We compute Equation ([26]) for two forms of the depen- 
dence of Fx on Mbh- For the first model, we assume 
a constant value of Fx = 2. For the second model, we 
assume that Tx depends on Mbh according to our best 
fit regression results, where Tx depends on Mbh ac- 
cording to Equation (Q2J) for M Bff < 3 x 10 8 M Q , and 
Tx depends on Mbh according to Equation (fT4|) for 
Mbh ^ 3 x 10 8 M Q . We ignore the intrinsic dispersion 
in r^f ■ In Figure [T4l we show the spectra computed from 
Equation (26]) for RQQs with M B h/M q = 10 7 , 10 s , 10 9 , 
and 10 10 . The dependence of the location of the peak in 
the disk emission on Mbh is clearly illustrated. 

In Figure [15] we compare the C( ox -Mbh regression re- 
sults with the dependence of a ox on Mbh expected from 
Equation ([26]) for both Tx~Mbh models. Under the thin 
disk approximation, a correlation is expected between 
a ox and Mbh, even if the fraction of bolometric luminos- 
ity emitted by the disk is independent of Mbh- However, 
our data are inconsistent with the assumption that fc- 
and Mbh are independent, given the thin disk approx- 
imation. Under the assumption that fjj and Mbh are 
independent, the oi ox -Mbh correlation is too flat, and a 
increase in the fraction of bolometric luminosity emitted 
by the disk with increasing Mbh is needed to match the 
steeper observed dependence of a ox on Mbh- Alterna- 
tively, if fo increases with increasing to, as we argue in 
§ 17.21 then a steeper a ^-MsH correlation also results if 
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Mbh and rh are correlated. In this case, if fn increases 
with increasing to, and if Mbh increases with increasing 
to, then fry will also increase with increasing Mbh, thus 
producing a steeper observed cx ox -Mbh correlation. 

To the extant that Equations (|23f (|26|) accurately ap- 
proximate the spectral shape of RQQs, our data imply 
that either the fraction of bolometric luminosity emitted 
by the disk increases with increasing Mbh, that Mbh 
and rh are correlated, or both. Some theoretical models 
have suggested that the fraction of bolometric luminosity 
emitted by the disk should depend on to, but b e relatively 
insen sitive to Mbh (e.g., ICzernv et al.1 120031 : iLiu et all 
2003). Therefore, while a significant dependence of fn 
on Mbh is not predicted by these disk/corona models, 
these models are still consistent with the interpretation 
that a M^/f-TO correlation is driving the steeper a ox - 
Mbh correlation. Unfortunately, without accurate esti- 
mates of to we are unable to distinguish between these 
two possibilities. 

A shift in the peak of the disk SED with Mbh may 
also explain the dependence of a ox on redshift observed 
by K07. K07 speculated that the observed hardening of 
a ox with increasing z at a given Luv may be due to a 
correlation between a ox and Mbh, manifested through a 
Mbh~z correlation. At a given Luv, an increase in Mbh 
will result in an increase in Lx relative to Ljjv, assum- 
ing that to is not strongly correlated with Mbh- This is 
because an increase in Mbh decreases the temperature 
of the disk, shifting the peak in the disk SED toward the 
red, and thus increasing the luminosity at 2500A. How- 
ever, since K07 investigated the dependence of a ox on z 
at a given Luv, the luminosity at 2500 A is held constant. 
Therefore, the overall disk emission must decrease in or- 
der to keep the luminosity at 2500A constant despite the 
increase in Mbh- As a result, an increase in Mbh at a 
given Luv will result in an increase in the X-ray lumi- 
nosity relative to the luminosity at 2500A. Because Mbh 
and z are correlated in our flux limited sample (e.g., see 
Figure [2]), an increase in z will probe RQQs with higher 
Mbh- As a result, RQQs will become more X-ray loud 
with increasing z, at a given 2500A luminosity. Con- 
sequently, deeper surveys that probe a greater range of 
Mbh should not see as strong of a correlation between 
Mbh and z, thereby reducing the magnitude of a a ox -z 
correlation. Indeed, investigations based on samples that 
span a greater range in luminosity do not find evidence 
for a, correlation betw een a ox and z (e.g., ISteffen et al.l 
l2006tlJust et al.ir2007f L qualitatively consistent with our 
interpretation of a a ox -z correlation. 

7.2. Dependence of a ox on rh 

We have found that a ox increases with increasing 
Ljjv I LEdd, and decreases with increasing Lx/LEdd- 
The mere existence of these correlations is not partic- 
ularly interesting, as we would expect that the ratio 
of optical/UV luminosity to X-ray luminosity would in- 
crease as the fraction of optical/UV luminosity rela- 
tive to Eddington increases, and vice versa for an in- 
crease in Lx/LEdd- However, the relative magnitude 
of these dependencies carries some information regard- 
ing the dependence of a ox on to. A correlation between 
a ox and Luv / LEdd implies that Luv/Lx increases as 
the quantity fuvfh increases, where fuv is the fraction 
of bolometric luminosity emitted at 2500A. Likewise, an 



anti-correlation between a ox and Lx / L Edd implies that 
Luv/Lx decreases as the quantity fxrh increases, where 
fx is the fraction of bolometric luminosity emitted at 2 
keV. If the fraction of the bolometric luminosity emit- 
ted by the disk, as compared to the corona, increases 
with increasing m, then we would expect a strong in- 
crease in Luv/Lx with the product fuv™, resulting 
from the dual dependency of Luv/Lx on fuv and rh. 
Furthermore, because the fraction of bolometric luminos- 
ity emitted by the disk should decrease with increasing 
fx, then, if the fraction of bolometric luminosity emit- 
ted by the disk increases with increasing to, we would 
expect a weaker dependence of Luv/Lx on the quan- 
tity fxfh- This is because an increase in to causes 
an increase in the disk emission relative to the corona 
emission, which will then work against the decrease in 
disk emission relative to the corona that results from 
an increase in fx- The end result is a weaker depen- 
dence of Luv/Lx on the product fxrh- Indeed, this is 
what we observe, where Luv/Lx °c (Luv / LEdd) 2 ' 5 and 
Luv/Lx oc {Lx/LEdd)^ 1 ' 5 ■ Therefore, we conclude that 
the disk emission relative to the corona emission increases 
with increasing rh. This is in agreement wi th some mod- 
els o f corona with a slab geometry (e.g.. ICzernv et all 
[19971 Uaniuk fc Czernvl 120001 : iMerloni fc Fabianl 120021 : 
ILiu et al.l l20031 where the a ox -rh correlation arises due 
to a dependency of the size of the corona on rh. 

Our result that a ox is correlated with Luv / LEdd and 
anti-correlated with Lx / LEdd is inconsistent with a con- 
stant bolometric correction to both the optical/UV and 
X-ray luminosities. Instead, an increase in Luv/Lx 
with increasing rh implies that the bolometric correc- 
tion depends on rh. Because we conclude that the frac- 
tion of bolometric luminosity emitted by the disk in- 
creases with increasing to, this implies that the bolo- 
metric correction to the optical/UV luminosity decreases 
with increasing to, while the bolometric correction to 
the X-ray luminosity increases with increasing rh. The 
direction of this trend is c onsistent with the results of 
IVasudevan fc Fabianl ((2007), who find that the bolomet- 
ric correction to the X-ray luminosity increases with in- 
creasing Lboi/ LEdd- Similarly, we have found evidence 
that the fraction of bolometric luminosity emitted by 
the disk depends on Mbh, therefore implying that the 
bolometric correction also depends on Mbh- Even if 
the fraction of bolometric luminosity emitted by the disk 
is independent of Mbh, the bolometric correction will 
still depend on Mbh because the location of the peak 
in the disk emission will shift toward longer wavelengths 
as Mbh increases. As Mbh varies, the luminosity at 
2500A probes a different region of the quasi-blackbody 
disk emission, thereby producing a dependence of bolo- 
metric correction on Mbh- 

7.3. Implications for Black Hole Feedback 

A significant amount of recent work suggests that ra- 
diative and mechanical feedback energy from AGN plays 
an important part in galaxy and s upermassive black hole 
coevolution (e.g.. iFabianl Il999t IWvithe fc Loebl 120031 : 
Hopk ins et al. I I2005TT Within the context of these mod- 
els, a nuclear inflow of gas, possibly the result of a galaxy 
merger, feeds the SMBH, thus igniting a quasar. The 
SMBH grows until feedback energy from the quasar is 
able to drive out the accreting gas, thus halting the ac- 
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cretion process. Hydrodynamic calculations of accretion 
flows have shown that the efficiency of the quasar in driv- 
ing an outflow depends on the fraction of energy emitted 
through he UV/disk co mponent as compared to the X- 
ray/corona component (Proga 2007). The disk compo- 
nent produces luminosity in the UV, which is responsi- 
ble for driving an outflow via radiation pressure on lines, 
whereas the corona component produces luminosity in 
the X-rays, which is responsible for drivi ng an o utflow 
via thermal expansion. Calculations by iProgal (|2007f ) 
have shown that radiation driving produces an outflow 
that carries more mass and energy than thermal driv- 
ing. If the efficiency of black hole feedback depends on 
the quasar SED, any dependence on Mbh and rh of the 
fraction of AGN energy emitted in the UV as compared 
to the X-ray has important consequences for models of 
black hole growth. 

Because we have found evidence that the fraction of 
bolometric luminosity emitted by the disk increases with 
increasing rh and Mbh, this implies that black holes be- 
comes more efficient at driving an outflow with increasing 
rh and Mbh- However, the olq^-Mbh correlation may 
be due to the combination of both a correlation between 
Mbh and rh, and a dependence of the location peak in 
the disk SED on Mbh- If the fraction of energy emitted 
by the disk only depends weakly on Mbh, as some the- 
oretical models have suggested (e.g., ICzernv et al.l 120031 : 
iLiu et al.l [2003h . the fraction of energy emitted in the 
UV will still decrease with increasing Mbh becuase the 
peak of the disk emission will shift away from the UV. 
In this case, at a given rh we would expect that black 
holes will become less efficient at driving an outflow with 
increasing Mbh- 

7.4. Dependence of Tx on Mbh and rh 

In this work we have also found evidence that Tx 
and Mbh, Ljjv / LEdd, and Lx/LEdd are not statisti- 
cally independent. Moreover, the dependence of Tx on 
black hole mass or Eddington ratio appears to follow 
a non-monotonic form, although the Tx~Mbh trend is 
weak compared to the dependency of Tx on Eddington 
ratio. For the Tx~Mbh relationship, the X-ray con- 
tinuum hardens with increasing black hole mass until 
Mbh ~ 3 x 10 8 M©, after which the X-ray continuum 
softens with increasing black hole mass. The opposite 
is true for the Tx-L uv /L E dd and T x -L x /L E dd trends, 
and further work is needed to confirm this result. Previ- 
ous studies have not seen this non-monotonic trend be- 
cause they have only employed the H/3 emission line, and 
therefore their samples have been dominated by low-z, 
\ow-Mbh sources. 

7.4.1. Selection Effects 

It is unlikely that the dependence of Tx on 
M B H,L uv /L E dd, or L x /L E dd is due to redshifting of 
the observable spectra range. If this were the case, then 
as Mbh increases, so docs z due to selection effects, and 
thus we would observe a decrease in Tx as the 'soft ex- 
cess' shifts out of the observed 0.3-7 keV spectral range, 
while the compton reflection component shifts into the 
observed spectral range. However, there are lines of 
evidence that suggest that the Tx correlation are not 
due to redshifting of the observable spectral region, and 
that at least some of the observed dependency of Tx on 



M B H,L uv /L E dd, and L x /L E dd is real. First, in § [O] 
we tested whether a sample of nine z > 1 test sources 
with Mbh <3x 10 8 Af© were better described by a re- 
gression fit using the z > 1.5, Mbh ^ 3 x 10 8 Af Q sources, 
or by a regression fit using the z < 1, Mbh ^3x 10 8 M© 
sources. We found that the test sources were better 
fit using the regression of similar Mbh, and therefore 
that the difference in the Tx~Mbh correlations primar- 
ily depends on Mbh- Second, similar trends at low 
redshift between Tx and Mbh or L^oi/LEdd have been 
seen in other studies that only analyze the hard X-ray 
spectral slope (typically 2- 12 keV, e.g., iPiconcelli et al.1 
2005; Shc mmer et al J 12006). and thus these studies are 
not effected by the soft excess. Third, the compton re- 
flection hump is unlikely to shift into the observabable 
spectral range until z ~ 1. However, the contribution 
to the inferred Tx from compton reflection at z > 1 is 
likely weak, if not negligible, as our z > 1 sources have 
Mbh ^ 10 8 Mq and are highly luminous, and therefore 
are expected to only have weak reflection components 
(iMineo et al.ll2000t iBallantvne et~aLll2001t iBianchi et all 
120071) . 

There are two scenarios in which the non-monotonic 
behavior of Tx with Mbh or Eddington ratio may be 
artificially caused by selection. We will focus on the Ed- 
dington ratio dependency, as it is the strongest; however, 
our argument also applies to Mbh- First, the intrin- 
sic dependency of Tx on Eddington ratio could be lin- 
ear with increasing intrinsic scatter at high L\, i/ LEdd- 
Then, an inferred non-monotonic trend would occur if we 
were to systematically miss quasars with high Lb i / L Edd 
and steep X-ray spectra. Alternatively, there could be no 
intrinsic dependency of Tx on Eddington ratio. In this 
case, we would infer a non-monotonic trend if we were 
to systematically miss quasars with steep X-ray spectra 
at low and high hboil i->Edd, and quasars with flat X-ray 
spectra at moderate L^oil 'LEdd- 

We do not consider it likely that the observed non- 
monotonic dependence of Tx on Eddington ratio is due 
solely to selection effects. K07 describes the sample se- 
lection for sources with Tx ■ With the exception of some 
of the z > 4 quasars, all sources from K07 were selected 
by cross-correlating the SDSS DR3 quasars with public 
Chandra observations. Almost all SDSS sources in K07 
had serendipitious Chandra observations, and therefore 
were selected without regard to their X-ray properties. 
K07 estimated Tx for all sources that were detected 
in X-rays at the level of 3a or higher. Therefore, the 
only additional criterion beyond the SDSS selection im- 
posed by K07 is the requirement that the source had to 
be detected in X-ray, which was fulfilled by 90% of the 
quasars; the undetected sources were slightly more likely 
to be found at lower redshift, probably due to the pres- 
ence of unidentified BAL quasars. As a result, the K07 
sample selection function is essentially equivalent to the 
SDSS quasar selection function. Because the SDSS se- 
lects quasars based on their optical colors, the most likely 
cause of selection effects is the optical color selection. 
There is evidence that Tx is correlated with the slope of 
the optical continuum, where the X-ray continuum flat- 
tens (hardens) as the optic al continuum steepens (soft- 
ens) (iGallagher et al.l[2005L K07) The SDSS selectio n 
probability is lower for red sources (jRichards et al.ll2006f) . 
so we might expect to systematically miss sources with 
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smaller Tx- However, for the two scenarios described 
above, this is opposite the trend needed to explain the 
^x-Lx/LEdd relationship, where we need to at least sys- 
tematically miss sources with larger Tx- Furthermore, 
the drop in SDSS selection efficie ncy with optical spec - 
tral slope only occurs at 2 < z < 4 (jRichards et alfe OOB) , 
thus we would expect a redshift dependence for this se- 
lection effect. As we have argued above, and in § 16.31 
the non-monotonic trends for Tx cannot be completely 
explained as the result of different redshift ranges being 
probed. 

7.4.2. Implications for Accretion Physics 

The dependence of Tx on Luv / LEdd and Lx/LEdd is 
likely due to a dependence of Tx on rh. If these Tx cor- 
relations were due to a dependence of Tx on fuv or fx , 
then we would expect opposite trends for Luv / LEdd and 
Lx/LEdd, as fuv and fx should be anti-correlated. The 
fact that the regression results for the T x~Luv I LEdd 
and T x~Lx / LEdd relationships are similar implies that 
T x depends on rh, and at most only weakly on fuv or 
fx. 

A non-monotonic dependence of Tx on rh is predicted 
from the accreting corona model of iJaniuk fe Czernvl 
(2000), as well as a non-monotonic d ependence of T x 
on the viscosity (|Bechtold et al.l [20031 ). In addition, Tx 
is expected to steepen with increasing optical depth (e.g., 
lHaardt fc Maraschilll99lL[l99llCzernv et al.ll200l . One 
could then speculate that the dependence of Tx on Mbh 
or rh is due to a non-monotonic dependence of the corona 
optical depth on rh, which may indicate a change in the 
structure of the disk/corona system at ~ 3 x 10 8 M Q 
or some critical rh. Recent work also suggests a non- 
monotonic de pendence of the optica l /UV spectra l slope , 
otuv, on rh (jBonning et all 120071 : iDavis et all l2007f ). 
From this work, it has been inferred that the optical/UV 
continuum becomes more red with increasing rh until 
Lboi/LEdd ~ 0.3, after which the optical/UV contin- 
uum becomes more blue with increasing rh. Assum- 
ing the bolometric corrections described in § 14.21 the 
turnover in the Tx -Lboi I L Edd relationsh ip also occurs at 
Lboi/LEdd ~ 0.3. iBonning et al.l (|2007f ) suggested that 
the turnover in spectral slope at Lboi/LEdd ~ 0.3 may 
be due to a change in accretion disk structure, where the 
inner part of the accretion di sk becomes thicker due t o 
increased radi a tion p ressure (|Abramowicz et aTl Il988f ) . 
IBonning et al.l (|2007t ) performed a simple approximation 
to this 'slim disk' solution and found that it is able to pro- 
duce a non-monotonic trend between optical color and 
Eddington ration. If the inner disk structure changes 
at high rh, this change could alter the corona structure, 
producing the observed trend between Tx and Edding- 
ton ratio. 

Unfortunately, current models for corona geometry 
make a number of simplifying assumptions, and do not 
yet predict a specific relationship between a ox , Tx, Mbh, 
and rh. Ide ally, full magneto- hydr o dynami c sim- 
ulations (e.g., |T2 g_\2hier^_et_aI| 120031: iTurnerl l200l 
iKrolik et alj 20051) that include accretion disk winds 
(e.g.. lMurrav et al Hl995l : lProga fc Kallmanll2004D should 
be used to interpret the results found in this work. 
However, MHD simulations have not advanced to the 
point where they predict the dependence of a ox and 
Tx on quasar fundamental parameters, but hopefully re- 



cent progress in anal ytical descriptions of the magneto- 
rotational instability (|Pessah et al.ll2006t I2007D will help 
to overcome some of the computational difficulties and 
facilitate further advancement. 

8. SUMMARY 

In this work we have investigated the dependence of 
a ox and Tx on black hole mass and Eddington ratio using 
a sample of 318 radio-quiet qu asars with X-r ay data from 
ROS AT (|Strateva et al.ll2005D and Chandra (jKellv et all 
2007), and optical data mostly from the SDSS; 153 
of these sources have estimates of Tx from Chandra. 
Our sample spans a broad range in black hole mass 
(10 6 < Mbh/Mq < 10 10 ), redshift (0 < z < 4.8), and 
luminosity (10 43 < Ai A (2500A)[erg s" 1 ] < 10 48 ). To 
date, this is the largest study of the dependence of RQQ 
X-ray parameters on M B h, Luv /L E dd, and L x /L E dd- 
Our main results are summarized as follows: 

• We show that a ox is correlated with Luv / LEdd 
and anti-correlated with Lx/LEdd- This result 
is inconsistent with a constant bolometric correc- 
tion being applicable to both the optical/UV lu- 
minosity and the X-ray luminosity. This result, 
when taken in comb i nation with recent work by 
IVasudevan fc Fabianl (|2007l ). implies that constant 
bolometric corrections can be considerably unreli- 
able and lead to biased results. Instead, we argue 
that Luv /Lx increases with increasing rh and in- 
creasing Mbh, therefore implying that the bolo- 
metric correction depends on rh and Mbh- 

• We performed a linear regression of a ox on 
luminosity, black hole mass, Luv / LEdd, 
and Lx/LEdd, and found significant evi- 
dence that a ox depends on all four quan- 
tities: Luv/Lx oc L^ ±om ,L uv /L x oc 



M 



0.43±0.06 



N2.57±0.45 



J UV 

L BH °, Luv/Lx OC {Luv I LEddY 

and Luv/Lx oc Lx/L^^^ ' 14 . The dependence 
of a ox on Luv may be due to the dual dependence 
of a ox on Mbh and rh. Because we have attempted 
to correct for the statistical uncertainties in a ox 
and the broad line estimates of Mbh, these results 
refer to the intrinsic relationships involving a ox 
and Mbh, and are not merely the relationships 
between a ox and the broad line mass estimates, 
M BL oc LiFWHM 2 . 

• A correlation between a ox and Mbh is expected 
from the fact that the peak in the disk emission will 
shift to longer wavelengths as Mbh increases, even 
if the fraction of the bolometric luminosity emitted 
by the disk does not change with Mbh- Using a 
simple model for RQQ spectra, we argue that the 
observed cx ox -Mbh correlation is steeper than that 
expected if both rh and the fraction of bolometric 
luminosity produced by the disk are independent of 
Mbh- The observed o: ox -Mbh relationship there- 
fore implies that either the fraction of bolometric 
luminosity emitted by the disk increases with in- 
creasing Mbh, that Mbh is correlated with rh, or 
both. 

• A correlation between a ox and rh is predicted from 
several models of 'slab'-type corona. We argue that 
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the weaker dependence of a ox on Lx / L Edd implies 
that Ljjv/Lx increases with increasing m. Consid- 
ering that the efficiency of quasar feedback energy 
in driving an outflow may depend on the ratio of 
UV to X-ray luminosity, a correlation between a ox 
and both Mbh and to has important consequences 
for models of black hole growth. In particular, if su- 
permassive black holes become more X-ray quiet at 
higher to, they will become more efficient at driving 
away their accreting gas, thus halting their growth. 

• Because of a possible nonlinear dependence of Tx 
on M B H,Luv/L E dd, or L x /L E dd, we performed 
seperate regressions for the black hole mass es- 
timates obtained from each emission line. We 
confirmed the significant dependence of Tx on 
Ljjy I Lftdd and Lx/LEdd seen in previous stud- 
ies as inferred from the broad line mass estimates 
based on the H/3 line; however, we also find ev- 
idence that the Tx correlations change direction 
when including the C IV line. In particular, for 
the H/3 sample, the X-ray continuum hardens with 
increasing Mbh, while for the C IV sample, the X- 
ray continuum softens with increasing Mbh- Sim- 
ilar but opposite trends are seen with respect to 
Luv/LEdd and Lx/LEdd, and we conclude that 
these relationships can be interpreted as resulting 
from a correlation between T x and to. Results ob- 
tained from the Mg II line were too uncertain to 
interpret. We analyzed two test samples to argue 
that this non-monotonic behavior is not due to the 
different rcdshifts probed by the two samples, or to 
problems with the estimates of Mbh derived from 
the two lines; the different trends may be due to the 
difference in Mbh probed by the two samples. A 
non-monotonic dependence of Tx on Mbh and/or 
to may imply a change in the disk/corona struc- 



ture, although a non-monotonic dependence of T x 
on to and the viscosity is predicted by some models 
of 'slab'-type coronal geometries. 
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Fig. 1. — The [Ljjy , Lx, z) distribution of our sample. Non-detections are denoted by red arrows. 
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Fig. 2. — The distributions of estimated Mbh, Ltd, and L}, i/ Lfidd, as a function of z for our sample. The data points with error bars 
in the left two plots are fictitious data points illustrating the typical error in Mg^ and L bo i/ L Edd . 
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Fig. 3. — The distribution of a x as a function of Ljjy / L,Edd (left) and Lx/LEdd (right)- The opposite dependence of ct x on 
Luv I LEdd and Lx/LEdd suggests that at least one of these quantities is not proportional to L^ a i/LEdd- As such, we do not employ 
bolometric corrections in this work, and instead compare directly with Luv / LEdd and Lx/LEdd- 
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Fig. 4. — Ratio of optical/UV to X-ray flux as a function of Luv , Mgjf , Ljjy / LEdd-, and Lx/LEdd- The solid lines denote the best fit, 
and the shaded regions contain 95% (2cr) of the probability on the regression line. The data points with error bars in the plots of a ox as a 
function of Meu-,Ljjv /LEdd, and Lx/LEdd are fictitious and illustrate the typical errors in each direction. 
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Fig. 5. — Posterior probability distributions for the absolute value of the correlation coefficients between a x and 
\ogL uv ,\ogM BH ,\ogL uv /L Edd , and \ogL x /L Edd (left) and for the intrinsic dispersion in a ox at a given L uv , M BH , L uv /L Edd , 
or Lx/L Edd (right). The thick solid line denotes the probability distribution for the two quantities with respect to Ljjy, the dashed line 
for the two quantities with respect to Mgfj, the dashed-dotted line for the two quantites with respect to Ljjy / L Edd , and the thin solid 
line for the two quantities with respect to Lx/L Edd . 




Fig. 6. — Distribution of the X-ray photon index as a function of estimated Mbh, Luv/LEdd, and Lx/L Edd . For clarity, error bars arc 
only shown on Fx, and we cut-off the one data point with estimated Fx > 4. While no obvious trends between Fx and Mbh> Ljjy / L Edd , 
or Lx/L Edd exist for the whole sample, there is evidence of opposite trends in Fx for the H/3 and C IV samples. 
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Fig. 7. — Distribution of the X-ray photon index as a function of estimated Mbh, Luv / LEdd, and Lx / LEdd for the individual emission 
lines. For clarity, error bars have been omitted, and we omit the one data point with estimated Fx > 4. While no obvious trends exist for 
the whole sample, there is evidence of opposite trends for the H/3 and C IV samples. 




Fig. 8. — Posterior probability distributions of the slope (top) and intrinsic dispersion (bottom) from a linear regression of Fx on 
log Mbh , log Luv /LEdd, and log Lx/L^dd- The solid lines mark the posterior for the regression using the H/3 sample, the dashed-dotted 
lines mark the posterior for the regression using the Mg II sample, and the dashed lines mark the posterior for the regression using the C 
IV sample. 
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Fig. 9. — Posterior distribution for the difference in slopes between the C IV and H/3 regressions of Fx on tubh- While there is no 
significant evidence that either the H/3 or C IV regression slope is different from zero, there is significant evidence that they are not the 
same, implying a nonmonotonic trend between Tx and Mbh- 
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between emission line, Mbh-, an d z , we have collected a sample of H/3 test sources (red asterisks) at high Mbh an d z, and a sample of C 
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IV sample (dashed line). The high-z H/3 sources are better described by the high-z C IV regression, and the low-z C IV sources are better 
described by the low-z H/3 regression, implying that the difference in slopes between the H/3 and C IV samples is not due to systematic 
differences in mass estimates derived from either line. 
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Fig. 12. — Same as Figure ITOl but for the lower Mbh an d higher z test sources from COSMOS. These test sources help break the 
degeneracy between the Mbh an d z present in our main sample. 
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Fig. 13. — Same as Figure fTTT but for a sample of high-z, low- Mb h test sources from COSMOS (red asterisks with error bars). The 
COSMOS sources are better described by the low-z, low-Mgu regression, implying that the difference in slopes between the H/3 and C IV 
samples is due to the difference in Mjg probed by the two samples, and not due to the redshift differences. 




Energy [eV] Energy [eV] 

Fig. 14. — Model RQQ spectra computed from Equations J23D — . assuming Fx = 2 (left) and a varying Fx with Mbh (right). The 
spectra are computed for a RQQ with Mbh/Mq = 10 7 (thick solid line), 10 8 (dashed line), 10 9 (dashed-dotted line), and 10 10 (thin solid 
line). In all cases we assume m = 0.2 and that /jj = 85% of the bolometric luminosity is emitted by the disk. The vertical lines mark 
the locations of 2500A and 2 keV. The dependence of the location of the peak in the disk emission on Mbh ls apparent, producing a 
correlation between cr ox and Mbh even if the fraction of bolometric luminosity emitted by the disk is independent of Mbh- 
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M BH / M SUN 

Fig. 15. — Dependence of o x on Mbh compu ted from Equations l!23|l - l(2"6"Tl, a ssum ing Fx = 2 (blue dashed line) and a varying Tx 
with Mbh (blue dot-dashed line). As with Figure IT4l we compute Equations d23D — J26ft assuming rh = 0.2 and fr> = 0.85. The predictions 
from the model RQQ spectra are compared with our observed data and the regression results, where the symbols and lines have the same 
meaning as in Figure [3] The Oo^-Mbh relationships predicted from assuming that fp is independent of Mbh are inconsistent with the 
observed ce ox -MBH relationship. 



TABLE 1 

Black Hole Parameters of the Sample 



a(J2000) <5(J2000) z logA/ SL /M \ogL x /L Edd * logL uv /L Edd h 

00 02 30.7 +00 49 59.0 1.352 9.20 ± 0.45 -2.43 ± 0.45 -1.44 ±0.45 

00 06 54.1 -00 15 33.4 1.725 9.08 ± 0.29 -2.58 ± 0.30 -1.07 ±0.29 

00 22 10.0 +00 16 29.3 0.574 7.96 ±0.32 -2.73 ± 0.32 -1.11 ±0.32 

00 27 52.4 +00 26 15.7 0.205 7.35 + 0.45 -2.42 ± 0.46 -1.26 ±0.49 

00 31 31.4 +00 34 20.2 1.735 9.17 + 0.29 -2.66 ± 0.29 -1.31 ±0.29 

00 50 06.3 -00 53 19.0 4.331 9.61 ± 0.41 -2.89 ± 0.42 -1.19 ±0.41 

00 57 01.1 +14 50 03.0 0.623 8.66 + 0.32 -3.79 ± 0.32 -1.66 ±0.32 

00 59 22.8 +00 03 01.0 4.178 9.21 + 0.38 -2.28 + 0.41 -0.74 + 0.38 

01 06 19.2 +00 48 22.0 4.437 9.18 + 0.38 -1.96 + 0.39 -0.46 + 0.38 
01 13 05.7 +15 31 46.5 0.576 9.39 + 0.32 -3.47 + 0.32 -2.46 + 0.32 
01 13 09.1 +15 35 53.6 1.806 9.19 + 0.29 -2.32 + 0.29 -1.13 + 0.29 
01 15 37.7 +00 20 28.7 1.275 9.25 + 0.45 -2.75 + 0.45 -1.81 + 0.45 
01 26 02.2 -00 19 24.1 1.765 8.96 + 0.29 -2.58 + 0.29 -1.01 + 0.29 



Note. — The complete version of this table is in the electronic edition of the Journal. The 
printed edition contains only a sample. 

a Logarithm of the ratio of vL v [2 kcV] to L^^d, where L^dd is calculated from the broad 
emission line estimate of Mbh, Mbl^° Logarithm of the ratio of vL v [2500A] to LEdd, where 
L*Edd is calculated from the broad emission line estimate of Mbh,Mbl- 
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TABLE 2 

Sources with H/3 and C IV Used for Testing the Tx-Mbh Correlations 



Quasar Name 


a 


5 


z 


Line logM SL L bo i/L Edd Opt. Ref. a 




X-ray Ref. b 




(J2000) 


(J2000) 




M 







PG 0026+129 


00 


29 13.7 


+13 16 03.8 


0.142 


C IV 


7.789 


1.768 


1 


1.96 c 


5 


Fairall 9 


01 


23 45.7 


-58 48 21.8 


0.046 


C IV 


7.760 


0.265 


1 


1.83 ± 0.06 


6 


PG 1202+281 


12 


04 42.2 


+27 54 12.0 


0.165 


C IV 


7.855 


0.359 


1 


1.76 + 0.07 


7 


PG 1211+143 


12 


14 17.7 


+14 03 12.3 


0.080 


C IV 


7.559 


1.380 


1 


2.06 + 0.05 


8 


PG 1247+267 


12 


50 05.7 


+26 31 07.7 


2.038 


H/3 


10.41 


0.379 


2 


2.23 + 0.10 


9 


Q1346-036 


13 


48 44.1 


-03 53 25.0 


2.370 


H/3 


9.946 


0.609 


3 


2.02 ± 0.17 


3 


MRK 478 


14 


42 07.5 


+35 26 22.9 


0.077 


C IV 


7.890 


0.498 


1 


2.41 + 0.07 


7 


PG 1630+377 


16 


32 01.1 


+37 37 50.0 


1.476 


H/3 


9.762 


0.569 


4 


2.20 ± 0.30 


10 


PG 1634+706 


16 


34 28.9 


+70 31 33.0 


1.334 


H/3 


10.27 


0.734 


4 


2.19 + 0.05 


10 


HE 2217-2818 


22 


20 06.8 


-28 03 23.9 


2.414 


H/3 


10.12 


0.807 


3 


1.97 + 0.06 


3 



References. - - (1) IKellv fc Bechtoldl 1120071) (2 ) IMcIntosh et alj U9991) (31 IShemmer et al.l (1200611 (4) INishihara et al.l |1997D 
(51 IQ'Neill et al.l 020091 (61 INandra et al.l (1199711 (71 IBrocksopp etali 1 1200611 (81 IReeves fc Turner! (1200011 Ol lPage et al.l 1120041 (101 
IPiconcelli et al.l 1200*511 

a Reference for the rest frame optical/UV data.* 1 Reference for the X-ray data. c [O'Neill ct al. (20051 do n °t report an error on Tx- 



TABLE 3 
Test Sources from COSMOS 



a 


8 


z 


Line 


log M BL 


Lbol/LEdd 


Px 


(J2000) 


(J2000) 






M Q 






09 58 48.8 


+02 34 42.3 


1.551 


C IV 


8.276 


0.131 


2.01 ± 0.11 


09 59 02.6 


+02 25 11.8 


1.105 


Mg II 


7.612 


0.025 


2.17 + 0.28 


09 59 49.4 


+02 01 41.1 


1.758 


C IV 


8.108 


0.719 


2.51 + 0.16 


10 00 50.0 


+02 05 00.0 


1.235 


Mg II 


7.692 


0.501 


2.50 + 0.13 


10 00 51.6 


+02 12 15.8 


1.829 


Mg II 


7.807 


0.131 


2.14 + 0.17 


10 00 58.9 


+01 53 59.5 


1.559 


C IV 


8.346 


0.172 


2.04 + 0.17 


10 02 19.6 


+01 55 36.9 


1.509 


C IV 


8.333 


0.177 


2.19 + 0.23 


10 02 34.4 


+01 50 11.5 


1.506 


C IV 


7.991 


0.941 


2.25 ± 0.12 


10 02 43.9 


+02 05 02.0 


1.234 


Mg II 


7.817 


0.303 


1.97 + 0.29 



